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Formation and structure of the microemulsion phase in two-dimensional ternary 

AB+A+B polymeric emulsions 

Dominik Diichs and Friederike Schmid 

Fakultat fur Physik, Universitdt Bielefeld, 
Universttatsstr. 25, 33615 Bielefeld, Germany 

We present an analysis of the structure of the fluctuation-induced microemulsion phase in a 
ternary blend of balanced AB diblock copolymers with equal amounts of A and B homopolymers. 
To this end, graphical analysis methods are employed to characterize two-dimensional configuration 
snapshots obtained with the recently introduced Field-Theoretic Monte Carlo (FTMC) method. We 
find that a microemulsion forms when the mean curvature diameter of the lamellar phase coincides 
roughly with the periodicity of the lamellar phase. Further, we provide evidence to the effect of a 
subclassification of the microemulsion into a genuine and a defect-driven region. 



I. INTRODUCTION 

Microemulsions arise when two immiscible substances 
are compatibilized by the use of a surfactant lowering the 
interfacial tension between the incompatible components 
dramatically. Originally introduced for mixtures of oil, 
water, and a surfactant, the term microemulsion applies 
equally to polymeric blends. Microemulsions appear in 
two forms: droplet and bicontinuous. Here we consider 
the bicontinuous case, which is observed in symmetric 
systems with balanced surfactants. In polymeric emul- 
sions, the compatibilizing copolymer is usually the most 
expensive ingredient. A good theoretical understanding 
of these system is therefore desirable from an applica- 
tion point of view, as well. Applications often depend 
crucially on the morphology of the blend, and bicontinu- 
ity is a particularly useful property, e. <?., for improved 
stiffness or conductivity. 

Traditionally, there have been three distinct coarse- 
grain approaches to the theoretical study of self- 
assembling amphiphilic systems^]: (a) Microscopic 
approaches [3, often built around lattice-gas or lattice- 
Boltzmann simulations p| B , IE El El Hi ; also by dissipative 
particle dynamics |gL llOl iflj or by standard Monte Carlo 
or molecular dynamics methods [Hi Il3. Il4 lla. IWL fl7L fl8L 
0. (b) Ginzburg-Landau models^UMlilllii EJ, 
leading to a much simpler description by means of only 
a few order parameters and coefficients which can be ob- 
tained from experiment. Polymeric blends have been 
described very successfully with a particular continu- 
ous density functional approach usually denoted self- 
consistent field theory (SCFT)|H |H The draw- 
back of this method is that it employs a mean-field ap- 
proximation and therefore neglects the very cause for the 
existence of microemulsions in these systems, (c) Mem- 
brane approaches for systems in which the solubility of 
the amphiphilcs is extremely low, such that they can form 
membranes by themselves |26|. This case is not considered 
here. 



As mentioned above, the addition of a compatibiliz- 
ing agent, here the AB copolymer, is essential to the 
formation of mesoscopically ordered phases in a pure A 
and B homopolymer mixture. The melts would oth- 
erwise phase-separate macroscopically below some rela- 
tively high temperature, which is typically above room 
temperature. Microemulsions, although by definition 
part of the disordered regime, do have a finite (meso- 
scopic) preferential length scale, and thus their existence 
depends on the presence of a compatibilizer, as well. 
Apart from the microscopic correlation length, £, which 
is the characteristic decay length of the correlation func- 
tion, a second length scale, q^ 1 , which describes the wave- 
length of oscillations in the correlation function, comes 
into play. The disorder line is now defined as the locus 
where q just vanishes. The onset of the microemulsion 
regime, on the other hand, is defined by the Lifshitz line, 
which is the locus where the peak in the structure factor 
is just found at nonzero wavevector q* . A microemulsion 
is thus disordered but not structureless [H Efij . 

In a recent paper, we introduced a new approach to in- 
corporating the effect of thermal fluctuations in field the- 
ories of polymer melts: the Field-Theoretic Monte Carlo 
(FTMC) method|3fJ. It is an extension of the earlier 
mentioned SCFT method. The system under investiga- 
tion was the ternary model system consisting of symmet- 
ric AB copolymers as well as A and B homopolymers, the 
homopolymers being a = 0.2 times as long as the copoly- 
mers. This system had previously been the ob ject of a se- 
ries of experiments conducted by Bates et al. |3lll32ll33| . 
Mean-field calculations of its phase diagram predict a Lif- 
shitz critical point where the disordered, lamellar, and 
phase-separated regions meet [3, • However, it could 
be demonstrated both in these experiments as well as in 
the FTMC simulations that the Lifshitz point is destroyed 
by thermal fluctuations and a channel of bicontinuous mi- 
croemulsion (B[iE) emerges in between the lamellar and 
phase-separated regions. Whereas the aim of that pa- 
per, Ref. l3fl was mainly to establish the validity of the 
FTMC method and to compare it with a related com- 
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plex Langevin method |36l 1371 138| . we here want to ex- 
amine more closely the formation and structure of the 
microemulsion phase. 

Neither the disorder line nor the total monomer Lif- 
shitz line correspond to any thermodynamic transition. 
Recently, Morkved et al. 's have found that neither 
line correlates well with the transition from a fully disor- 
dered mixture at higher temperatures to a well-developed 
B/iE. Nevertheless, they found that dynamic light scat- 
tering provided a clear signal for this transition, which 
they proposed to be the homopolymer/homopolymer Lif- 
shitz line originally introduced by Holyst and Schick 1401 . 
In a similar experimental system, Schwahn et al. |4l| 
found evidence for the existence of three different regimes 
in the disordered phase, a "disordered blend", a "mi- 
croemulsion" , and a "disordered copolymer" region. The 
three regimes differ from each other by the value of the 
peak wavevector q* in the structure factor - zero in the 
disordered blend, large in the disordered copolymer, in- 
termediate in the microemulsion. On the theoretical side, 
Holyst and Przybylski^^ performed Monte Carlo simula- 
tions of a Ginzburg-Landau model for the lamellar phase 
in copolymers and showed that topological fluctuations 
change the monotonic decay of the off-specular scattering 
intensity, introducing the topological Lifshitz line. One 
might suspect that such fluctuations affect the character- 
istic lamellar distance to higher values. Experiments |43| 
as well as simulations |44j indeed reveal that the peak 
wavevector q* is ~ 15 — 20% lower at the order-disorder 
transition than predicted by SCFT. However, this effect 
can also be explained by local chain stretching (3 E3- 

With our approach, we can separate the effect of fluc- 
tuations from that of local chain rearrangements. It 
combines the advantages of SCFT, which accounts in 
full for the chain connectivity but neglects fluctuations, 
with those of Ginzburg-Landau simulations, which in- 
clude fluctuations but make approximations for the in- 
trachain correlations. Our results are qualitatively sim- 
ilar to those of Morkved et al. and Schwahn et al. 
Our main result in this paper is the subclassification of 
the microemulsion region into a disordered and a genuine 
regime. For the transition from the lamellar to this dis- 
ordered microemulsion phase we do not find a markedly 
lowered q*. In the genuine regime, q* differs from the 
mean field value. Indeed, this is our criterion for the dis- 
tinction of the two. It should be interesting to examine 
to what extent, if any, our subclassification corresponds 
to Morkved et al. 's [3^ distinction of poor and good 
microemulsions, or to Schwahn et al. 's distinction of 
disordered copolymers and microemulsions. 

In the present work, we have used parallel Cray archi- 
tectures to simulate AB+A+B melts on two-dimensional 
lattices of size 48x48. The high demands on computing 
power have so far limited us to two dimensions. It will be 



straightforward to adapt the method to three dimensions 
once this becomes feasible from a computational point 
of view. This paper is organized as follows. In Section 
II, we shall briefly outline the theoretical model under- 
lying FTMC. For a more detailed description, we refer 
the interested reader to our earlier paper[3(|. Section III 
contains the results and discussion. In this context, our 
evidence is presented to the effect of a subclassification 
of the microemulsion. We conclude in Section IV with a 
summary. 



II. THE FIELD-THEORETIC MODEL 



In this section, we briefly present the model underly- 
ing FTMC for the ternary system under consideration. 
For details, see Ref. |3(1 We study a mixture of ha 
homopolymers of type A, ub homopolymers of type B, 
and uab symmetric block copolymers in a volume V. 
The polymerization index of the copolymer is denoted by 
AT, and the corresponding quantities for the homopoly- 
mers are denoted as by Na = Nb — aN. We con- 
sider the case of a symmetric copolymer, wherein the 
fraction of A monomers in the copolymer is / = 1/2. 
We restrict our attention to the concentration isopleth, 
where the homopolymers have equal volume fractions 
4>ha = 4>hb = 4>h/2. The monomeric volumes of both 
A and B segments are assumed to be identically equal to 
l/po- On this model we impose an incompressibility con- 
straint. We model the effective interactions between seg- 
ments as Flory-Huggins local contact interactions and use 
a Gaussian chain model, which corresponds to perfectly 
flexible polymers. With these assumptions, the canonical 
partition function of the system can be written as 

Z c oc I DW- I DW + exp[-H c (W-,W + )] (1) 



with 

H C (W-,W+) = C 



X N 



drWt 



drW+ 



(2) 
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Po Rd 



2a 



(3) 



In this paper, all lengths are expressed in units of the un- 
perturbed radius of gyration, R g = b(N/(2d)) 1 ^ 2 , where 
d is the space dimension. The parameter C in the above 
equations, which occurs as a global prefactor to He acts 
as a Ginzburg parameter such that in the limit C — oo 
the partition function is reduced to its saddle point 
and the mean- field solution becomes exact. In Eq. 
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Qa, Qb, and Qab denote the single chain partition func- 
tions for the A, B, and AB chains, respectively, in the 
potential fields W- (r) and W+ (r) . Note that W— is con- 
jugate to the difference in A and B densities, m, and W+ 
to the total density, </>. Moreover, W- is real, whereas 
W+ is imaginary thereby rendering Hq complex. We 
here employ a partial saddle point approximation in W+ , 
which reduces He to real values. 

The single chain partition functions can be expressed 
in terms of the Feynman-Kac formulae as: 



Qi 



dv qi(r, z/j), 



(4) 



where the propagators % satisfy diffusion equations. 
From these propagators, we can calculate density opera- 
tors, 4>a and ipB, from 

Mr) = — / dsq A B(r,s)q AB (r,l-s) (5) 

Qab Jo 



V4>h 
2aQ A 



pCt 

/ ds q A (r,s)q A {r,a- s) 
Jo 



and a similar equation for the (\>b- The densities which 
correspond to the experimentally measurable quantities 
are the averages over 4>a,b, *■ e., 4>a.b — (4>a,b)- For 
the calculations in the present work, however, we have 
treated 4>a and 4>b as instantaneous densities. We then 
calculated time averages over parameters calculated on 
them. Strictly speaking, the quantities 4>A, b are visual- 
izations of the spatial distributions of fields. The time 
averages so obtained should nevertheless reflect the es- 
sential structural properties of the system. 

For a detailed description of the FTMC method used 
to obtain the configuration series that serve as the input 
data in this work, see Ref. |3(]. 



III. RESULTS AND DISCUSSION 

Fig. ^ shows the fluctuation-corrected phase diagram 
of the ternary melt, as obtained in Ref. jsjj, at dimension- 
less polymer number density C = 50, which is an inter- 
mediate value and has been used throughout this work. 
Compared with the corresponding mean-field diagram, 
one discerns (a) a shift in the order-disorder transition to 
higher segregation strengths, x^V, (b) a slight shift in the 
same direction of the transition between the disordered 
and phase-separated regions, and (c) the emergence of a 
cusp-like region of a microemulsion phase. Note that this 
diagram was calculated in two dimensions, due to numer- 
ical constraints. Likewise, the simulations carried out in 
the present work were in two dimensions only. While 
a full three-dimensional analysis will likely yield a more 
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Homopolymer Lifshitz line 
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FIG. 1: Phase diagram of the ternary A+B+AB blend (a = 
0.2). The solid lines show the mean-field phase diagram, which 
features a three-phase (L+A+B) coexistence region reaching 
up to the Lifshitz point. The circles give locations of fluctu- 
ation corrected phase boundaries at C = 50 (from Ref. |3(|). 
The thick dotted lines are guides for the eye. The order- 
disorder transition is weakly first-order, but the coexistence 
region is too small to be displayed. The thin lines in the dis- 
ordered region indicated the total monomer Lifshitz line (dot- 
ted), the homopolymer- homopolymer Lifshitz line (dashed), 
and the line of (g£) 2 = 3 (dashed-dotted), as obtained from 
mean-field theory. 



narrow microemulsion "channel" than that shown in Fig. 
^ the good qualitative agreement with experiment indi- 
cates that essential aspects of the phase behavior are cap- 
tured in two dimensions already. Thus it is justified to 
proceed with further analysis. 

In addition to these results taken from our earlier pa- 
per, Ref. l30l Fig. n also displays various Lifshitz lines 
calculated within the mean-field approximation. The to- 
tal monomer Lifshitz line is found at the homopolymer 
volume fraction 461 



1 



l + 2a 2 



(0) 



and is independent of the incompatibility parameter xA. 
The homopolymer-homopolymer Lifshitz line, on the 
other hand, is independent of a, and determined by the 
equation 



XN = 



$ff(l - 



(7) 



It was obtained by generalizing the calculation of Holyst 
and Schick[40| for homopolymer/copolymer length ratio 
a = 1 to arbitrary a. Note that in contrast to the case 
a = 1 studied by Holyst and Schick, the homopolymer- 
homopolymer Lifshitz line is rather close to the total 
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monomer Lifshitz line at a = 0.2. This is because short 
homopolymers swell the copolymer blocks, whereas the 
longer homopolymers in systems with a = 1 are expelled 
from the copolymer rich regions. 

Finally, Fig. also shows the line (q(,) 2 = 3, 

which distinguishes between microemulsions with weak 
and strong ordering tendency in confined geometries|47|. 
Here £ is the correlation length and q the wavevector q of 
oscillations in the correlation function. The line (q£,) 2 = 3 
was calculated by expanding the inverse mean-field struc- 
ture factor S^q) -1 up to fourth order in q, i. e., approx- 
imating S(q) by the Teubner-Strey form|48| . 



S(q) 



1 



(8) 



+ gq* + cq* 

and then determining the value of %iV where g 2 /log — 

We start by examining the structure factor of our two- 
dimensional melts. It can be calculated from the Fourier 
transform of the density correlation function, which in 
turn can be obtained from the fluctuating W A (r) = 
W+ (r) + W- (r) field via a formula derived in Ref . 
Thus, the structure factor 5(q) is the Fourier transform 
of 



1 

V 
with 



- / dr [(p(r )p(v + r)> - (p(r )) (p(r + r)>] , (9) 



2C* 



(p(r)p(r'))-(p(r))(p(r')) = — 6(r-r>) 



(10) 



4C 



[{W A {v)W A {v')) - (W A (r)) (W A {v'))\ 



(XN) 2 

In Fig. |3 structure factors for lamellar phases 
and strongly and weakly structured microemulsions are 
shown. The lamellar structure factor ( Fig. [21 a) features 
the hallmark double peak of a striped pattern. Fig. 
b) illustrates that in a microemulsion, anisotropy is lost 
yet a preferential length scale does exist, as evidenced by 
a ring-shaped region of maxima. As we progress deeper 
into the disordered phase at lower \N, the ring becomes 
less pronounced, until it is almost indistinguishable from 
a a shapeless low-amplitude noise ( Fig. [21 c). 

The above procedure to calculate the structure factor 
of a configuration is quite expensive from a computa- 
tional standpoint as it involves integrations over the en- 
tire course of a simulation. However, if we are to analyze 
geometrical patterns, it is not necessary to use the full- 
fledged structure factor. We can see the characteristic 
features of a configuration from the Fourier transform of 
its density distribution, as well. In this work, we have 
therefore used time averages of quantities derived from 
the Fourier transform of <p A in analyzing morphological 
properties. 



(a) 



(c) 





I S (q) | 
400 




FIG. 2: Structure factors at C = 50 (a) Lamellar phase at 
4>h = 0.7 and \N = 12.5. (b) Strongly structured microemul- 
sion at <f)H = 0.82 and xN = 12.5. (c) Weakly structured 
microemulsion at <f>H = 0.82 and \N = 10. 
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As a starting point, we need to find the preferential 
length scale of a density distribution (f>A (r) . This is equiv- 
alent to finding the preferential wave vector of its two- 
dimensional Fourier transform, F(q). To this end, we use 
the square averaged over all angles, defining 



Fo(q) 



1 

2^ 



271 



#|F(q)| 5 



(11) 



which is now one-dimensional. 



(a) 



JCN=12.5 
%N=12 
XN = 11.5 
XN=11 
XN= 10.5 
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(b) 




FIG. 3: (a) Equilibrated averages of Fo(q) at C = 50, 
<J>h = 0.82, and different values of \N ■ The arrow indicates 
the mean-field location q\j F of the maximum in the structure 
factor in the disordered phase, (b) Teubner-Strey fits to Fo(q) 
for selected values of \N- The inset shows snapshots of <j>A 
for the same \N values (12.5, 11, and 8 from left to right). A 
lattice point is painted black ("A") if 0.51 < (j>A < 1, grey if 
0.49 < <j) A < 0.51, and and white ("B") if < 4> A < 0.0.49. 

From Fig. in which Fo(q) is displayed for var- 
ious \N at ()>h — 0.82, it can be seen that F (q) has 
a pronounced peak for %iV ~ 11. Below \N ~ 11 
the curves become markedly flatter at low q. Note that 
the mean-field transition between the disordered and the 
lamellar phase occurs at \N = 11.4 (for = 0.82). 
Above this point, the position of the (hardly discernible) 
maximum of Fo(q) agrees well with its mean- field value, 



q*MpR g = 1.203 for = 0.82. In the more structured 
region of higher xiV, the peak moves to lower values of 
q, reflecting the fact that the corresponding mean-field 
length scale, the lamellar distance, increases with increas- 
ing \N. Fig. |21 b) shows fits of selected curves to the 
Teubner-Strey form ( Eq. ©) ^H- The fits are reason- 
able at high and low values of x_/V, and less satisfactory 
in the vicinity of the mean-field transition. 

The technicalities of calculating Fo(q) require q to be 
binned in a histogram and thus approximated by discrete 
values (in practice, by multiples of 2ir/L, where L is the 
box length). In addition, Fig. [3] illustrates that our data 
for Fo(q) are strongly scattered, mainly due to the statis- 
tical error. This makes the q coordinate of the maximum 
of Fo(q) a poor candidate to be taken as the preferential 
wave vector. A smoother procedure is to use 



fdq q F (q) 
JdqF Q (q) ■ 



(12) 



If Fo(q) has a pronounced well-defined maximum at 
nonzero q* ( e. <?., \N > H-5 in Fig. q roughly 

coincides with this maximum. Otherwise, it slightly over- 
estimates q* ( e. g., qoR g — 1.3 at %A = 11., qoR g = 1.6 
at \N — 8). Now, the preferential length scale, L , is 
defined as 



T 

-^o •— — • 

q 



(13) 



Further, we define the mean curvature diameter, Dc, 
of the boundaries of A and B microdomains in a black- 
and-white image similar to those in the inset of Fig. 
01b), in which a lattice point is painted black ("A") if 
0.5 < (f> A < 1, and white ("B") if < 4> A < 0.5: 



Dn := 2 





[ ds 


dt 


2" 






L C J 




ds 





(14) 



L c is the sum of all contour lengths of the microdomain 
boundaries, and t is the tangent vector at a given co- 
ordinate s along the contour. For a detailed explana- 
tion of the algorithm used to calculate Dc, see Ap- 
pendix A. Note that the snapshots shown here also con- 
tain grey pixels, corresponding to a balance density, i.e., 
0.49 <4> A < 0.51. 

In Figs. 01 El configurational snapshots 49] and 
plots of equilibrated averages of Lq and Dc are displayed 
for various homopolymer volume fractions, (f>n, at %A = 
12.5, 12, and 11.5. The deviation of Lq from the corre- 
sponding mean-field values, Lmf, is shown in the inset. 
The Monte Carlo simulations were run up to approxi- 
mately 1.5 million Monte Carlo steps, including equili- 
bration times of up to a few 100,000 steps. A Monte 
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FIG. 4: Characteristic length scales at x N = 12.5, C = 50. 
(a) Snapshots at 4> H = 0.74, 0.78, 0.8, and 0.82 (first row from 
left), and <j> H = 0.84, 0.86, 0.88, and 0.9 (second row from 
left), (b) Preferential length scale, Lo (circles), and curvature 
diameter, Dc (squares), in units of R g , vs. &h, averaged 
over the equilibrated parts of simulations on a 48 x 48 lattice. 
The solid line shows the mean- field lamellar distance Lmf for 
comparison. Filled and empty symbols correspond disordered 
and lamellar initial conditions, respectively. The inset shows 
the difference of Lo and Lmf- 



Carlo step includes one attempted random local incre- 
ment of W- per pixel, within ranges that were chosen 
such that the Metropolis acceptance rate was 35 %. As a 
general rule, the more disordered the configurations, the 
shorter was the equilibration phase. 

In all three cases, i.e., \N — 12.5, 12, and 11.5, Dc is in 
good approximation constant over an extended 4>h range 
( Figs. 01 E3 and . This is plausible because the curva- 
ture is induced by the copolymers, which are located pre- 
dominantly at the interfaces between microdomains. It 
should thus be independent of the amount of homopoly- 
mers. Nevertheless, Dc gets somewhat smaller as X-ZV 
is decreased; the higher temperatures that correspond to 
lower xN facilitate the bending of the microdomain in- 
terfaces. 

At xiV = 12.5, coming from the lamellar phase at low 
4>h , the curvature radius Dc becomes comparable in size 
to the preferential length scale, Lq, as well as the mean- 
field length scale, Lmf, around <f>u rs 0.75, which is when 




FIG. 5: Same as Fig. H for X N = 12. As in Fig. H 
the snapshots in (a) correspond to <j>H = 0.74, 0.78, 0.8, and 
0.82 (first row from left), and <j> H = 0.84, 0.86, 0.88, and 0.9 
(second row from left). 



the lamellae begin to break up. This is in good quan- 
titative agreement with the location of the fluctuation- 
corrected order-disordered transition of Fig. ^ At \N 
= 12 and 11.5, Dc is well below both the mean-field peri- 
odicities, Lmf, an d the preferential length scales, L , in 
the <j>jj range examined: these configurations are in the 
disordered phase. We have thus identified the mechanism 
by which fluctuations generate the microemulsion. 

The preferential length scale Lq and the mean-field 
length scale Lmf are identical within error bars at small 
homopolymer concentrations $h- Upon increasing 
at \N = 12.5 and 12, one notes that Lo begins to devi- 
ate from Lmf at a certain . This signals the departure 
from defect-driven behavior and the onset of a more gen- 
uine morphology within the microemulsion phase. The 
effect is stronger for \N = 12.5 than for \N = 12. At 
\N = 11.5, it is barely noticeable anymore, and Lq coin- 
cides within error bars with Lmf- However, the higher 
X-ZV, the more smeared out Lq becomes. 

As pointed out already, a microemulsion is character- 
ized by the existence of a preferential length scale within 
the disordered phase 0, |4g. In this sense, all configu- 
rations displayed in Figs. 0] and are microemul- 
sions. The previous result, however, suggests a more di- 
versified classification into (a) a defect-driven region and 



(a) 




(b) 




FIG. 6: Same as Fig. gjfor x N = U.5 The snapshots in (a) 
correspond to (f>u = 0.82, 0.83, 0.84, and 0.85 (first row from 
left), and cj} H = 0.86, 0.87, 0.88, and 0.89 (second row from 
left). 



(b) a "genuine" microemulsion. We shall denote them 
D/iE and G/iE, respectively. In the first region, the mi- 
croemulsion retains the characteristic length scale of a 
fluctuation-free lamellar phase, and the main effect of 
the fluctuations is to introduce lamellar defects, which 
destroy the lamellar order. Typical configurations are 
the first snapshots (top row) in Figs. 0] and [5] a). 
In the second region, the interfaces themselves fluctuate, 
as can be observed in the images in the bottom row of 
Figs. 01 and a). This causes an effective increase of 
the characteristic length scale (or decrease in the charac- 
teristic wavevector). Deeper in the disordered phase, at 
low \N, the distinction between lamellar defects and in- 
terfacial fluctuations is less obvious, one rather observes 
general structured, but unspecific disorder. In view of the 
fact that the characteristic wavevector coincides with the 
mean-field periodicity or, in the disordered phase, q* MF ( 
Fig. |3J), we classify it as D[iE. 



IV. SUMMARY AND CONCLUSION 



In this paper, we have investigated the structure of 
the disordered microemulsion phase in two-dimensional 
symmetrical ternary AB+A+B polymeric blends in dif- 
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FIG. 7: Substructure of the microemulsion. DfiE: defect- 
driven, GfiE: genuine microemulsion morphology, 20: phase- 
separated region. The solid lines show the fluctuation- 
corrected phase diagram of Fig. the dotted lines the 
mean-field phase diagram, and the dashed line is just a guide 
for the eye. 



ferent regions of the phase diagram. Not surprisingly, 
the structure of the microemulsion is strongest at param- 
eters X-/V and <pH where the mean-field approximation 
would predict an ordered phase, i. e., where the dis- 
order is brought about solely by fluctuations. We have 
corroborated our earlier result |30| regarding the mech- 
anism that underlies the formation of the microemul- 
sion phase in ternary AB+A+B polymeric blends: the 
lamellae break up when the curvature diameter of the 
microdomain boundaries becomes comparable to the pe- 
riodicity of the lamellar phase |5pj. We have further shown 
that the preferential length scale in the system deviates 
from its mean-field equivalent in a parameter subspace 
of the microemulsion only, demonstrating that the mi- 
croemulsion region is divided into a defect-driven and a 
"genuine" part. In the resulting redrawing ( Fig. 0) of 
the phase diagram Fig. ^ we have indicated the ap- 
proximate location of the conjectured transition between 
the two, as marked by the occurrence of nonzero values 
for Lo — Lmf m Fig- H]andE] Comparing Fig. [7]with 
the various mean-field Lifshitz lines indicated in Fig. ^ 
we find that both the total monomer Lifshitz line and the 
homopolymcr-homopolymer Lifshitz line are far from the 
region of interest here. The numerical results presented 
in Fig. |3 confirm that this is still true for the true 
(fluctuation-corrected) total monomer Lifshitz line. We 
have not investigated the homopolymer-homopolymer 
structure factor. Since the short homopolymers in our 
system (a = 0.2) swell the copolymers, the homopolymer- 
homopolymer Lifshitz line is not very far from the total 
monomer Lifshitz line. 

It is possible that the G[iE phase coincides with the 



"microemulsion" (jaE) phase indicated in the experimen- 
tal phase diagram of Schwahn et al. |4lJ . Strictly speak- 
ing, however, the two systems cannot be compared di- 
rectly to each other, because our simulations were car- 
ried out in two dimensions only. We hope that three- 
dimensional calculations will become feasible in the fu- 
ture. 

We thank M. Matsen, V. Ganesan, and G. Fredrickson 
for fruitful discussions. This work was supported by the 
Deutsche Forschungsgemeinschaft (Germany). The sim- 
ulations were carried out on the CRAY T3E of the NIC 
institute in Jlilich. 



APPENDIX A: GRAPHICAL ANALYSIS 



2+ \'2 1+2 V2 



V2 1 2 V2 2 2 V2 

SHE 




FIG. 8: Local contributions to the circumference. The current 
(black) pixel from step (3) here is painted grey. All diagrams 
are modulo 7r/2 rotations. 



In order to analyze configuration snapshots, it is often 
advantageous to take the 4>a distribution and convert its 
continuous values (between and 1) to black-and-white 
bitmaps with white pixels for < <pA(x,y) < 0.5, and 
black pixels for 0.5 < 4>a{x,u) < 1. From these images, 
we can then extract the following parameters. 



1. Algorithm to calculate the circumference, L c 

(1) Besides the original bitmap, Bi, define another 
bitmap, i?2, and paint it white, i.e., set all pixels to 
0. 

(2) Copy all border pixels from B\ to B^. A border 
pixel is defined as one which is black (i.e., has a 
value of 1) and whose four nearest neighbors have 
at least one white pixel among them. 

(3) Scan through B2 from the upper left to the lower 
right corner and for each pixel check whether it is 
black. If so, proceed as follows: 

(3a) If any pair of adjacent nearest-neighbor (NN) 
and next-nearest-neighbor (NNN) pixels of the 
current pixel are both black, paint the NN 
pixel white. This is done to prevent ambigui- 
ties in the subsequent steps. In practice, how- 
ever, this case is quite rare and barely changes 
the end result. 

(3b) Check the distribution of black pixels among 
the nearest and next-nearest neighbors accord- 
ing to Fig. |S] and add the specified numbers 
to the circumference. 

(4) Once the entire lattice has been scanned, divide the 
result by 2 to account for double counting in step 
(3b). 



2. Algorithm to calculate the curvature diameter, 

D c 

The curvature diameter of the (combined) circumfer- 
ence, L c , of all black areas is defined as 





f ds 


dt 


2" 


L c J 




ds 





with t the local (normalized) tangent vector. To calculate 
the quantity 



[ ds 


dt 






ds" 



(A2) 



we pursue the following strategy: 

(1.1) Make a border pixel bitmap as above. 

(1.2) Scan the bitmap and stop when a black pixel is 
found. 

(1.3) Start a data array, (xpUj), for the coordinates 
(numbered j > 0) of a new line, i > 0. 

(1.4) Look for a (nearest or next-nearest) neighbor of the 
current pixel. In most cases, there will be only one 
possibility to continue the line. Otherwise, an arbi- 
trary choice is made. 

(1.5) Remove the current pixel from the bitmap. Make 
the neighboring pixel from (4) the current pixel. 

(1.6) Repeat (1.4-5) until no neighbor is found for the 
current pixel. Then terminate that line. 

(1.7) Repeat (1.2-6) until the bitmap has been cleared. 
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(II. 1) For each line i, = (xj,yj), as determined in 

(1.1-7) that is longer than 3 points, set j = 1. 



(II. 2) Calculate |dt/ds| 2 according to 



dt 
ds" 



dt x = 



dty = 



2 = dt|+dt| 

ds 2 ' 



'3+2 - Xj 





x j) 2 + (Vi+2 - Vj) 2 




XjJ r \ Xj — i 




- Xj-i) 2 + (Vj+i -Vj-i) 2 








Xj) 2 + (Vj+2 - Vj) 2 




Vj+i - Vj-i 



(A3) 
(A4) 

(A5) 



y/(x j+1 - Xj^Y + (y j+1 - Vj^f ' 



dx = X 3+2 - + X] ~ Xj-l 



(A6) 



ds = \J dx 2 + Ay 2 . 



(A8) 



(11. 3) Increase j. Repeat (II. 2-3) until the end of the line 
is reached. 

(11.4) Repeat (II.1-3) for all lines. 

(11.5) Apply 



(A9) 



Note that this algorithm produces a slightly different 
(smaller) value for L c than the one presented in the last 
section. This because the length of each line is effectively 
truncated by three points. L c here acts only as a nor- 
malizing factor and for consistency should be calculated 
as indicated. If interested in L c itself, one should use the 
algorithm of the previous section. 



dy = 



y ]+ 2 - yj+i + yj - Vj-i 



(A7) 
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